NASA Contractor Report CR-187127 


/4 /£5 7 ^^ 

p. no 


Investigation, of Advanced Counterrotation 
Blade Configuration Concepts for High Speed 
Turboprop Systems 


Task IV - Advanced Fan Section Aerodynamic Analysis 
Computer Program User’s Manual 


Andrew J. Crook and Robert A. Delaney 
Allison Gas Turbine Division of General Motors 
Indianapolis, Indiana 


November 1992 

NASA 

Prepared for 

Lewis Research Center 

Under Contract NAS3-25270 


(N AS A-CR- 187127) INVESTIGATION OF N93-18702 

AOVANCEO COUNTERROTATION BLADE 
CONFIGURATION CONCEPTS FOR HIGH 

SPEED TURBOPROP SYSTEMS. TASK 4t Unclas 

ADVANCED FAN SECTION AERODYNAMIC 
ANALYSIS COMPUTER PROGRAM USER'S 

MANUAL Report, Jul. 1991 - Jul . Hl/07 0145784 

1992 (General Motors Corp. ) 70 p 




Preface 


This manual was prepared by Andrew J. Crook and Robert A. Delaney of the 
Allison Gas Turbine Division, General Motors Corporation, Indianapolis, IN. The 
work was performed under NASA Contract NAS3-25270 from July, 1991 to July, 1992. 
The modification of the flow analysis for multiple splitter fan/engine geometries was 
performed by Andrew J. Crook. The Allison program manager for this contract was 
Robert A Delaney. The NASA program manager for this contract was Christopher 
J. Miller. 


Acknowledgements 

The authors would like to express their appreciation to the following personnel 
who contributed to this program: 

Dr. Edward J. Hall for his suggestions concerning the solver and the many helpful 
contributions to this document, 

Dr. Christopher J. Miller for his review of the program, 

UNIX is a trademark of AT&T 

IRIS is a trademark of Silicon Graphics, Inc. 




PREC*0»NG page BLANK' NOT FILMED 



TABLE OF CONTENTS 


1. SUMMARY - 1 

2. INTRODUCTION 3 

3. ADPAC-APES PROGRAM INITIALIZATION AND DEMON- 
STRATION 7 

3.1 Introduction 7 

3.2 Extracting the Source Files . . . 7 

3.3 Compiling the Source Code 8 

3.4 Running the Distribution Demonstration Test Case 10 

4. ADPAC-APES PROGRAM GENERAL DESCRIPTION ...... 15 

4.1 Parameter /Array Sizes 13 

4.2 Input File Description - Basic Input Section 16 

4.3 Input File Description - Blade Row File Name Section 24 

4.4 Input File Description - Extra Input Parameter Section 26 

4.5 Output File Description 31 

4.6 Subroutine Description 35 

5. ADPAC-APES SOLUTION TECHNIQUES 41 


IV 



5.1 Single/Multiple Blade Row Procedure 41 

5.2 Solution Startup Techniques 49 

REFERENCES 53 

APPENDIX A. FULLPLOT : PostScript X-Y Plotting Routine 55 

APPENDIX B. ADPAC-APES DISTRIBUTION LIST 65 


V 


1HTSNV. 


UWMA-* 


ftE «• 



1. SUMMARY 


The purpose of this study is the development of a three-dimensional Euler/Navier- 
Stokes flow analysis for fan section/engine geometries containing multiple blade rows 
and multiple spanwise flow splitters. An existing procedure developed by Dr J. J. 
Adamczyk and associates [1] at the NASA Lewis Research Center was modified to 
accept multiple spanwise splitter geometries and simulate engine core conditions. The 
procedure was also modified to allow coarse parallelization of the solution algorithm. 
This document is a user’s manual for the code, developed under Task IV of NASA 
Contract NAS3-25270. 

The numerical solution is based upon a finite volume technique with a four stage 
Runge-Kutta time marching procedure. Numerical dissipation is used to gain solu- 
tion stability but is reduced in viscous dominated flow regions. Local time stepping 
and implicit residual smoothing are used to increase the rate of convergence. Mul- 
tiple blade row solutions are based upon the average-passage system of equations. 
The numerical solutions are performed on an H-type grid system, with meshes being 
generated by the system (TIGG3D) developed earlier under this contract. The grid 
generation scheme meets the average-passage requirement of maintaining a common 
axisymmetric mesh for each blade row grid. 
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The analysis was run on several geometry configurations ranging from one to five 
blade rows and from one to four radial flow splitters. Pure internal flow solutions were 
obtained as well as solutions with flow about the cowl/nacelle and various engine core 
flow conditions. The efficiency of the solution procedure was shown to be the same 
as the original analysis. / _ j 
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2. INTRODUCTION 


This document contains the Computer Program User’s Manual for the ADPAC- 
APES (Advanced Ducted Propfan Analysis Codes- Average Passage Engine Simula- 
tion) program developed by the Allison Gas Turbine Division of the General Motors 
Corporation under Task IV of NASA Contract NAS3-25270. The objective of this 
task is development of a three-dimensional flow analysis tool for advanced fan section 
and turbofan engine geometries such as the NASA/GE Energy Efficient Engine seen 
in Fig. 2.1 . The tool is able to compute steady flow solutions about geometries with 
any number of blade rows and any number of axisymmetric radial flow splitters. The 
tool simultaneously computes the flow through every flow path both internal and ex- 
ternal, which includes the flow through the fan, and about the fan cowl and nacelle, 
from upstream to downstream of the engine. When the domain is extended in this 
manner, engine performance can be determined entirely by the analysis tool. Effects 
of engine core flow can also be simulated. 

User instructions for setting up the computer program and running it for a 
demonstration geometry are found in Chapter three of this document. Chapter four 
covers some details on the subroutines, input and output of the program. Instructions 
and advice for running the code is given in Chapter five. 

This flow analysis tool was developed from a code entitled VSTAGE which was 
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developed by John J. Adamczyk of the NASA Lewis Research Center [1] . The 
user is referred to the documentation for that code for additional information on the 
operation of the code. 
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Figure 2.1: NASA/GE Energy Efficient Engine cross section 
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3. ADPAC -APES PROGRAM INITIALIZATION AND 
DEMONSTRATION 


3.1 Introduction 

This section describes the commands necessary to extract the ADPAC-APES 
source code from the standard distribution and run a complete test case for a flow 
solution about a ducted fan with a part span shroud. The standard APES distribution 
is normally a compressed tar file which can be decoded into the various parts by a 
sequence of commands on any standard UNIX system. The sequence listed below is 
intended to guide the user through the setup from the standard distribution and the 
running of the demonstration case. 


3.2 Extracting the Source Files 


The APES programs are distributed as a compressed tar file named 

apes. tar. Z 

The first step necessary to extract the ADPAC-APES programs is to uncompress the 
tar file with the command: 

uncompress apes.tar.Z 
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This operation replaces the compressed file apes.tar.Z with an uncompressed file 

apes. tar. 

The next step is to extract the individual files and directories from the apes. tar 
file. This process will create a subdirectory named apes in the current directory, so 
the user should move the apes. tar file to a suitable initial directory. Once the tar 
file is properly placed, the ADPAC-APES distribution may be extracted with the 
command 

tar xvof apes. tar 

(On some systems tar xvf apes. tar may be sufficient.) Execution of the command 
Is -1 will verify that the APES directory has been created. 

3.3 Compiling the Source Code 

After extracting the source files, the user is naturally interested in compiling the 
source files for execution. A UNIX-compatible Make facility is provided for the APES 
program. The Makefile which governs the compilation process is necessarily machine- 
dependent and requires that the user select from one of a number of preconfigured 
systems. The systems which are immediately available are: 

iris Silicon Graphics Iris workstation 
Cray Cray Computer Inc. supercomputer 
aix IBM Aix operating system UNIX workstation 
If no system is specified, then the iris system is assumed. The machine dependence 
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of the compilation process is inherently tied to the use of the Scientific DataBase 
LIBrary routines for binary file input/output. 

In order to begin the compilation, it is first necessary to enter the APES directory 
with the command: 

cd apes 

At this point, several files and directories will be available. By entering the command 
Is, a listing of the individual directories can be obtained. The output of the Is 
command will look something like: 

demo/ manual/ report/ sdblib/ src/ 

A description of each of these listings is given below: 

demo This directory contains all the files necessary for generating a demon- 
stration solution with the ADPAC-APES code. 

manual This directory contains the LaTeX source code for this manual. If 
LaTeX is installed on your system, it is possible to reproduce this 
document (excluding figures) with the command latex manual. The 
resulting device independent file manual. dvi may then be converted to 
PostScript or previewed on screen through a number of widely available 
routines. 

report This directory contains the LaTeX source code for the final report out- 
lining the technical details of the ADPAC-APES codes. If LaTeX is 
installed on your system, it is possible to reproduce the final report (ex- 
cluding figures) with the command latex report. Again the resulting 
file report. dvi may converted to a print output or previewed. 
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src This directory contains all the FORTRAN source code for the APES 
program. 

sdblib This directory contains the various machine-dependent files for the 

Scientific DataBase Library routines [2]. j 

In order to compile the APES code, enter the src directory with the command: j 

cd src 

The compilation is begun by issuing the command { 

make system j 

i 

* 

1 

where system indicates the current computing platform (iris, cray, or aix) described 
above. 

3.4 Running the Distribution Demonstration Test Case 

Once the make facility has properly completed compiling the APES source code, 
it is possible to run the test case provided. It is recommended that the test case be 
run to verify proper compilation and extraction of the APES distribution. 

In order to run the demonstration case, it is necessary to begin in the demo 
directory. At this point, the demo directory may be entered by issuing the command 

cd ../demo 

In the demo directory, an Is command will reveal: 

demo . input 1 gridl . demo mover* 

demo.input2 grid2.demo script. flip* 
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Free Stream Boundary 



Figure 3.1: Demonstration case conditions and geometry 

The demonstration case is a two blade row fan section with a cowl and a free stream 
outer boundary. Figure 3.1 shows the basic geometry and flow situation. The file 
demo. input 1 is the input file for the first blade row (rotor) of the demonstration case, 
and demo.input2 is an input file for blade row number two (stator). The blade row 
specific grids for the rotor and the stator are files named gridl.demo and grid2.demo , 
respectively. The mover* and script. flip* files will be explained later. To start the 
demonstration case, the user can issue the command 

../src/apes < demo.inputl > outputl 

This command executes the APES program, directing the input to be read from 
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the file demo.inputl and directing the output to the file outputl. When the program 

is finished, the Is command will reveal several new files: 

> [ 

I 

axil .new bfl.new flowl. new f lowl .new .hist 

f lowl .new .maxhist f lowl .new . p3dabs f lowl .new .p3drel f lowl .new .rmshist 

outputl | 

The file axil. new is the axisymmetric solution, the file bfl.new contains the blade 
body forces, and the file flowl .new contains the flowfield restart for the first blade 
row. All three files are needed to continue a multiple blade row calculation, while 
only flowl. new is needed for a single blade row calculation. A file containing the 
residual history data for the current run is flow 1 .new. hist. The residual history data 
files contain the maximum density time derivative (flowl. new. maxhist) and the root 
mean square average of the density derivative (flowl .new. rmshist) for each iteration. 

A PostScript plot of the history data can be generated using the FULLPLOT tool 
provided in the distribution (see the Appendix for a description of FULLPLOT). 

Before the plot file can be generated, the FULLPLOT code must be compiled by 
issuing the command: 

f77 -o ../src/fullplot .. /src/fullplot. f 

When this has been done, a plot file of the residual history is produced by issuing the 
command: 

../src/fullplot < flowl. new. hist 

The resulting file fort. 15 can be previewed or plotted using a PostScript compatible 
printer. The files flowl .new. maxhisi and flowl. new. rmshist contain the maximum 
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density residual and the average residual, respectively. The purpose of these files will 
be explained later. Finally, the files flowl.new.pSdabs and flowl.new.p3drel contain 
the absolute and relative PL0T3D flowfield data interpolated to the grid nodes, 
respectively. These files are formatted as PLOT3D multiple grid (mgrid=l) binary 
files and can be used in combination with the grid gridl.demo to examine the solution 
with PLOT 3 D. 

To initialize and run the second blade row the user can issue the command: 

../src/apes < demo.input2 > output2 

The input file directs the code to initialize the stator blade solution with the ax- 
isymmetric solution of the rotor. The rotor blade’s body forces (read in from the 
file bfl.new ) are incorporated into the flowfield and the stator specific solution com- 
mences. When the program is finished, several more files will appear in the directory: 

axi2.new bf2.new flow2.new flow2 .new .hist 

f low2 .new .maxhist f low2 .new ,p3dabs f low2 .new.rmshist output2 

These files contain the restart data, residual history data, and PL 0T3D flowfield data 
for the stator blade row solution. 

At this point one cycle through all the blade rows in this case (one flip ) has 
been completed. For this demonstration case, several flips are needed before a steady 
solution is obtained. See [3] or the Solution Techniques section of this manual for 
more information of flips. To set up for another flip, the user issues the command. 

mover 
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This executes the file mover * which contains several file manipulation commands 
which are simply executed in order. The commands rename the flowfield restart, 
axisymmetric solution, and body force files. Other commands append the blade row 
residual history files to a general residual history file so that a continuous record of 
the solution residuals can be plotted. [ 

To execute a complete flip, the input files to APES must be changed so that f 

the solution calculation can continue from its restart files. When this is done the 
commands: j 

../src/apes < demo.inputl >> outputl 

../src/apes < demo.input2 >> output2 j 

1 

mover 

can run the demonstration case through one flip. All of these changes and commands 
have been included in the file script. flip*. This file executes the program with the 
correct input files for each blade row and then rearranges all the necessary files to set i 

up for the next flip. Thus, one flip can be accomplished by issuing the command: 

script.flip 

After several flips (about 10) the solution should be converged. The solution is judged 
as converged when the root mean square average of the density residual (i.e., the 
quantity in the file withe the suffix new.rmshist) has been reduced by a factor of 
10~ 3 from its initial value. For more information about flips and convergence can be 
found in the Solution Techniques Section. 
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4. ADPAC-APES PROGRAM GENERAL DESCRIPTION 


4.1 Parameter/ Array Sizes 

It is necessary to estimate the maximum array sizes and other limit parameters 
prior to the compilation process. Most of the arrays in the routines of the code 
are dimensioned according to the PARAMETER statements in the file parameter. inc 
which contains the line: 

PARAMETER ( IMX=100, JMX=50, KMX=50, MSP=6 ) 

These parameters dimension numerous arrays found in the subroutines. The 
PARAMETER variable are defined as. 

IMX Number of streamwise grid lines + 1 
JMX Number of spanwise grid lines + 1 
KMX Number of circumferential grid lines + 1 
MSP Maximum number of spanwise flow splitters 
Before running the code, the user should know the size of the grid and make sure 
the array sizes are adequate. The PARAMETER variables can be larger than defi- 
nitions above but errors and peculiar results will occur if the arrays are dimensioned 
below the required parameter. 
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4.2 Input File Description - Basic Input Section 


A sample input data file for ADPAC-APES is shown in Fig. 4.1. The beginning 
of the input file consists of an arbitrary number of header lines (for user information) 
which is terminated by a line beginning with a character. Following this, title 
information is read in and will be printed out in the output file. There can be up to 
nine lines of title information. If more than nine lines are present, then only the first 
eight lines of information and the last line will be output from the code. As with the 
header section, the title information is terminated by a line beginning with a 
character. Beyond this line, the file follows the structured format which is seen in 
Fig. 4.1, beginning with the 4l H — MACH” characters. The data is structured as to 
the order it is read in but there is no field oriented format in which the input must 
be organized. The field markers seen in Fig. 4.1 are only to assist the user as to the 
order of the input data. Most input data is read in as a real number except where 
indicated. A brief description of the variables used in the data file is given below. 
Some of the variables in the input file are not used in the current code but their place 
is kept so that continuity with other flow solvers is maintained. Values are needed for 
these inactive variables since they hold places and are read in before active variables. 
The inactive variables are italicized in the descriptions below. 

VARIABLE DESCRIPTION 


MACH Freestream Mach number. 

GAMMA Specific heat ratio (cp/cy). 

PEXIT Exit static pressure. This is applied at the hub and 
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User information about this particular blade row or running condition 


+ TITLE ( up to 9 lines ) 

nrMnwq^RA^ION INPUT FILE FOR A GENERIC FAN ROTOR 

this'fan rotor is one blade row in a three blade row calculation 

THERE IS ONE SPLITTER (CORE/BYPASS ) AND A CORE CONDITION 
REPRESENTING THE TURBINE EXIT IS SIMULATED 
THE GRID IS 150 BY 50 BY 50 FOR THIS BLADE ROW 
TAKE OFF CONDITIONS ARE MODELED FOR THE SYSEM 

+ --MACH — ++-GAMMA-++-PEXIT-++-OMEGA— ++--ADVR--++--DOCT--++-BCFAC-- + 

0 3500 1.4000 0.0000 -4500.000 0.00 1.0 1.0000 

+ __CFL ++— VIS2— ++--VIS4— ++— QFIL— ++ 3C ++-BLDROW-++- DTAXI--+ 

_ 4 0000 2 00 1.000 0.0 0.0000 1.0 1.00 

+ -FNCMAX — H+--REST — + + -"~SAVE--+ + ”FISTEP-++-FNPRNT-+- , --ROWMAX-++-FPRAXI- + 

100 0 -1.0 -1.0 1.0 0.0 3.0 00.0 

+-REFINI-++-REFINJ-++-REFINK-f +--EPSX ++ EPSY ++ EPSZ--++-FPRBF--+ 

, 0 1.0 1.0 2.00 2.00 2.00 0.0 

+ -RFINII- + + -RFINJJ-++-RFINKK-++— FAXI — ++-FAXIMX-++-FR3DMX-++— FTURB- + 

RFIN l o i o 1.0 0.0 99999.0 0.0 0020.0 

+-FSTRT3-++-FEXIT3-++ — FSET — ++-FINVIS-++-JBASE ++ JTIP--++-f wf c j--+ +_ fWi.ck--+ 

0 c 0.0 0.00 0.0 1.0 5 j . 0 1.0 1.0 

+ __ D t AM „+ + --TREF-+ + -PREF-++-RGAS-++-PR— ■ -++-PRT— ++- trans-++-itrans-+ 
5 5750 518.7 2116.8 1716.48 0.72 0.90 0.0 1.0 

+ -OMEGS1-++--FCKL H FCHT--+-* OMEGS2-++“"FCTL--T+--FCTT--+-i--f iwf s--++ _ f iwf e--+ 

-4500.00 9.0 150.0 0.0 0.0 ^ °‘ C £ ^i° 

low. rest .fori blade_force . fori axi_f row . *.orx g. id. for- 

axi_f low . f or 2 
axi_f low . f or3 
axi_ flow . new! 

-LSPLIT-+ — LINLET-+ 

1 1 


lew. rest .for2 blade_f orce . for2 

low. rest ,for3 blade_f orce . f or3 

low. rest .newl blade_f orce . newl 

+-NZXITS-+ — LCORE--+ — NMAVG — + 

1 1 2 


grid . f or2 
grid. for3 


PEXIT2-+ — PEXIT3-+ ( EXIT PRESSURES ) 

ENGINE CORE SIMULATION VARS 


f-PEXITi-+ 

1.310 

+ -ICORIN-+— ICORSX-+— JCORE—+— FCORE— + 

141 141 19 1.0 

+ -PSCCRE PTCORE- + — TTCORE-^- — ANCORE-+ ( AN 


IN DEGREES ) 


1.0 1.0 

- — I A VG — t J 1 AV — + J 2 A V 

123 1 19 

137 21 59 

1--NSPLIT- + 

1 

+--JSPL — + — ISPLE--+--ISPTE — + 
19 75 141 

+-LPTFLG-+— LTTFLG-+ — LANFLG-+- 
0 0 0 

+ — PTIN — + TTIN — + ANIN — +■ 

2116.0 520.0 10.0 

2326.0 420.0 20.0 

2116.0 320.0 30.0 


0.0 


( MASS AVG . PLANE STREAMWISE S SPANWISE BOUNDS ) 

( NUMBER OF RADIAL SPLITTERS ) 

( SPLITTER SPANWISE & STREAMWISE BOUNDS } 

-NINPTS-+ ( INLET PROFILE PARAMETERS ) 

-3 

— RPOS — + 

0 . 10 
0.50 
0 . 90 


Figure 4.1: Sample input data file for ADPAC-APES 
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OMEGA 

ADVR 

DUCT 


BCFAC 

CFL 


integrated radially outward for shrouded geometries 
(DUCT=1.0), or applied at the outer exit boundary and 
integrated radially inward for propellers (DUCT=0.0), 
to satisfy radial equilibrium. 

(this parameter may not be used when domain has multiple 
exits, see exit pressure definition in Extra Input Section) 

Rotational speed (revolutions per minute) of the 
particular blade row. Positive value indicates rotation 
counter-clockwise as one is looking down the machine axis. 

Not active in this version. 

Internal flow duct option: 

if = 0.0, external flow options are utilized (unshrouded blade with free 
stream outer boundary), 

if = 1.0, an internal flow is assumed (compressor). 

Solid wall boundary conditions are applied at the outer 
boundary rather than the characteristic far field 
condition. 

Not active in this version. 

Time step parameter: 
if < 0 then this is the CFL number used 
to determine the calculation time step for 
local time stepping (i.e., for each cell), 
if > 0.0 then this is the CFL number used 
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to determine the calculation time step without local 

time stepping (i.e., most restrictive cell CFL sets time step throughout 

domain). 

VIS2 Second order damping coefficient (« 1.0 — 3.0) 

(divided by 4 in the code) 

VIS4 Fourth order damping coefficient (« 1.0 — 2.0) 

(divided by 64 in the code) 

QFIL Artificial dissipation trigger: 

if < 0.0 then the standard dissipation scheme is used, 
if > 1.0 then the dissipation is directionally scaled. 

(this eigenvalue scaling is sometimes usefull for highly stretched grids) 
BC Exit boundary condition trigger: 

if = 0.0, characteristic boundary condition based on Pi 0 i is used, 
if = 1.0, characteristic boundary condition based on mass flow is used. 

BLDROW Blade row 7 parameter. This value determines which blade row 7 
is being calculated during a multiple blade row 7 solution. 

For single blade row calculations, this should be = 1.0. 

DTAXI Not active in this version. 

FNCMAX Maximum number of time steps to be performed. 

REST Restart option parameter: 

if = —1.0 a restart file and body force file are read, 

(multiple blade row configuration) 

if = 0.0 no restart file, initialize variables in code, 
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SAVE 


FISTEP 

FNPRNT 


ROWMAX 


FPRAXI 


if = 1.0 a restart file is read. 

(single blade row configuration) 

Save restart file option parameter: 

if = -1.0 restart files and body source terms are output 
at the end of a run (multiple blade row configuration), 
if = 0.0 no restart file is output at the end of the run, 
if = 1.0 a restart file is output at the end of the run. 

(single blade row configuration) 

Number of iterations performed between recalculation of the time step. 

Flowfield output parameter: 

if = 0.0 no flowfield data printed out, 

if = 1.0 inlet flow (to leading edge) quantities written to output file, 
if = 2.0 flowfield written out in PLOT3D multiple 
grid, binary Q file format. 

Maximum number of blade rows in the current configuration. 

For single rotation, = 1.0. For multiple blade rows, 
this parameter indicates that additional grids and 
body source terms are to be read. 

Print axisymmetric solution, 

if = 0.0 axisymmetric flow quantities are not printed, 

if = 1.0 axisymmetric flow quantities (ASCII) are written to the output 

file. 


REFINI Not active in this version. 



REFINJ 

REFINK 

EPSX 

EPSY 

EPSZ 

FPRBF 


RFINII 

RFINJJ 

RFINKK 

FAXI 

FAXIMX 

FR3DMX 

FTURB 

FSTRT3 

FEXIT3 

FSET 


Not active in this version. 

Not active in this version. 

Implicit residual smoothing coefficient in the streamwise 
direction (« 2.0 is a typical value). 

Implicit residual smoothing coefficient in the spanwise 
direction (« 2.0 is a typical value). 

Implicit residual smoothing coefficient in the circumferential 
direction (% 2.0 is a typical value). 

Print trigger for body forces: 

if = 0.0, body forces are not printed, 

if = 1.0, body forces (ASCII) are written to the output file. 

Not active in this version. 

Not active in this version. 

Not active in this version. 

Number of axisymmetric solutions, not active in this version (0.0). 
Number of 3D passes before axisymmetric solutions, not active in this 
version (99999.0). 

Number of reduced domain 3D passes, not active in this version (0.0). 
Number of iterations performed before recalculating 
the turbulent viscosity. 

Not active in this version (0.0). 

Not active in this version (0.0). 

Flowfield restart parameter: 
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FINVIS 


JBASE 

JTIP 


FWFCJ 


FWFCK 


if = 0.0, standard restart procedure used, 

if = 1.0, restart using the previous axisymmetric solution, 

if = 2.0, restart using a piecewise construction of the 

flowfield from the axisymmetric solution of the other blade rows, 

if = 3.0, restart with a redefined uniform flowfield from 

inlet to the blade row leading edge and from the trailing 

edge to the exit. 

Viscous/inviscid solution trigger. This variable determines whether a 
viscous or inviscid solution is performed: 
if = 1.0, inviscid; 
if = 0.0, viscous. 

User defined base of the blade row. 

User defined tip of the blade row. This overrides the 
spanwise tip index found in the grid file. 

NOTE: The parameters JBASE and JTIP specify 
the bounds from which flow passes over and under 
the blade (a periodic surface), regardless of whether or 
not the mesh defines a periodic surface at those bounds. 

Wall function trigger for the endwalls. 

if = 0.0, shear stress near endwalls found from velocity gradient, 
if = 1.0, wall function used to determine shear stress 
near the endwalls and splitters. 

Wall function trigger for blade surfaces. 
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DIAM 

TREF 

PREF 

RGAS 

PR 

PRT 

TRANS 


ITRANS 


OMEGSl 

FCHL 

FCHT 


if = 0.0, shear stress near blades found from velocity gradient, 
if = 1.0, wall function used to determine shear stress 
near the blade surfaces. 

The rotor diameter or the value of the nondimensionalizing 
factor used in the input grid (feet). 

Reference total temperature (degrees Rankine). 

Reference total pressure (pounds per square foot). 

Gas constant (foot-pounds per slug degree Rankine). 

Gas Prandtl number. 

Turbulent Prandtl number (0.9 recommended). 

Fraction of the chord from the leading edge that transition 
on the suction surface begins. 

NOTE: All other surfaces have turbulent boundary layers. 

The suction surface indicator for the transition model: 

if = 1.0, the suction surface is the first circumferential index of the 

grid (usually if OMEGA is negative), 

if = 2.0, the suction surface is the last circumferential index of the grid 
(usually if OMEGA is positive). 

Rotational speed of a user specified section of the hub (RPM). 

Grid index for the hub or spinner leading edge (see Fig. 4.2) rotating 
at OMEGSl (used only in viscous calculations). 

Grid index for the hub trailing edge (see Fig. 4.2) rotating at OMEGSl 
(used only in viscous calculations). 
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OMEGS2 

FCTL 

FCTT 

FIWFS 

FIWFE 


Rotational speed of a second user specified section of the hub (RPM). 
Grid index for the hub leading edge (see Fig. 4.2) rotating at 0MEGS2 
(used only in viscous calculations). 

Grid index for the hub trailing edge (see Fig. 4.2) rotating at 0MEGS2 
(used only in viscous calculations). 

Axial grid index at which the wall function model starts. 

Axial grid index at which the wall function model ends. 


4.3 Input File Description - Blade Row File Name Section 


The section following the basic input parameters contains the names of files 
needed to restart a solution and the names of files written out following the solution 
iteration procedure. That section of the input file with suggested file names is repro- 
duced below. For each blade row in the domain (= BLDROW parameter) a line in 
the input file is required. This line contains four names of files pertinent to that blade 
row solution. The first file name is the flowfield restart for that blade row. The second 
name is for the blade row body force file. The third name is for the axrsymmetric 
flowfield, and the fourth is for the grid particular to that blade row. 


flow. rest. fori 
f low. rest. for2 
flow. rest .for3 
flow. rest .newl 


blade_f orce . fori 
blade.f orce . f or2 
blade_f orce . f or3 
blade_f orce .newl 


axi_f low. fori 
axi_flow.for2 
axi_flow.for3 
axi_f low. newl 


grid. fori 
grid.for2 
grid.for3 


Following the line(s) pertaining to blade row restart data, one line is required 
for the outputted restart data. Names for the updated restart flowfield, blade forces, 
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FCHL FCHT FCTL FCTT 


Figure 4.2: Fan section displaying hub rotation regions/parameters 
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and axisymmetric solution for the current blade row are on this line. The names can 
be up to thirty characters long but all four (or three) names must not take up more 
than one eighty column line. The names should be separated by blanks or commas. 


4.4 Input File Description - Extra Input Parameter Section 

- • ■ - • •• - * ' Sat : " - 

The third and last section of the input file contains parameters important to mul- 
tiple splitter calculations, engine simulations, and some user specified post processing. 
The parameters are read in a free format and are both integer and real. Explanations 
of the parameters follow and a number of these are illustrated graphically in Fig. 4.3. 

VARIABLE DESCRIPTION 1;" • 


NEXITS Number of flowfield domain exits (int). If one or more flow splitters 
extends to the exit, the domain exit is broken up (Fig. 4.3) and bound- 
ary conditions must be independently specified for each exit section. 
NEXITS specifies how many non-core (see LCORE below) exit sections 
there are. 

LCORE Core boundary condition flag (int). 

If = 0, no core is simulated. 

If = 1, a core inlet (and possibly exit) is simulated. A core inlet is 
treated as a flowfield domain exit with appropriate boundary condi- 
tions. A core inlet is sometimes defined at the last axial grid index (see 
Fig. 4.3) instead of a standard exit because the output file includes core 
mass flow in its iteration history when a core is specified. 
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Parameters: NEXITS = 2 NMAVG = 1 



ISPTE (1st and 2nd splitters) 


Parameters: NEXITS = 1 LCORE = 1 NMAVG = 0 NSPLIT = 2 



Figure 4.3: Computational domain highlighting some input parameters 
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NMAVG The number of mass averaged streamwise computational planes (int). 

This is a post-processing parameter to specify the number of stream- 
wise grid planes (I index=const) for circumferential mass averaged and 
global mass averaged data. 

LSPLIT Flow splitter definition flag (int). 

If = 0, no splitters are defined by the user/ 

If = 1, splitter bounds are input by the user. 

LINLET Inlet condtion definition flag (int). 

If = 0, standard inlet conditions are used. 

If = 1, user defines inlet conditions/profile. 

PEXITl Exit static pressure for the first exit section (real). This is applied 
at the lower endwall and integrated radially outward, or applied at 
the outer exit boundary and integrated radially inward (if the exit 
section has a free stream outer bound), to satisfy radial equilibrium. 
See PSCORE definition. 

PEXIT2 Exit static pressure for the second exit section etc. 

ICORIN Streamwise grid index for the core inlet (int). 

ICOREX Streamwise grid index for the core exit (int). 

JCORE Spanwise grid index for the top of the core passage (int). The bottom 
of the core passage is the first spanwise index. 

FCORE Not currently used (real). 

PSCORE Core inlet static pressure at the hub (real). 

Note: If a core is specified then PSCORE sets that domain exit pressure 
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PTCORE 


TTCORE 


ANCORE 


IAVG 

JlAV 

J2AV 

NSPLIT 

JSPL 

ISPLE 

ISPTE 

LPTFLG 


LTTFLG 


and the next (spanwise) domain exit pressure will be set by PEXIT or 
PEXIT1. 

Core exit total pressure (real). 

This parameter has no effect on the flowfield if ICOREX is greater 
than the last streamwise grid index. 

Core exit total temperature (real). 

This parameter has no effect on the flowfield if ICOREX is greater 
than the last streamwise grid index. 

Core exit flow angle in degrees from axial (real). 

This parameter has no effect on the flowfield if ICOREX is greater 
than the last streamwise grid index. 

Streamwise grid index for a mass average data plane (int). 

Spanwise grid index where the mass average plane begins (int). 
Spanwise grid index where the mass average plane ends (int). 

Number of flow splitters (int). 

Spanwise grid index at the lower surface of the splitter (int). 
Streamwise grid index at the leading edge of the splitter (int). 
Streamwise grid index at the trailing edge of the splitter (int). 

Inlet total pressure trigger (int): 

If = 0, a constant total pressure (reference pressure) profile is used. 

If = 1, a spanwise total pressure profile is input (lbs/square foot). 

If = 2, a spanwise total pressure profile is input (nondimensional Pt). 
Inlet total temperature trigger (int): 
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If = 0, a constant total temperature (reference temp) profile is used. 
If = 1, a spanwise total temperature profile is input (degrees Rankine). 
If = 2, a spanwise total temperature profile is input (nondimensional 


LANFLG 


NINPTS 


PTIN 

TTIN 

ANIN 

RPOS 


Tt). 

Inlet absolute swirl angle trigger (int): 

If = 0, a constant swirl angle (0 deg from axial) profile is used. 

If = 1, a spanwise swirl angle profile is input (deg from axial). 

If = 2, a spanwise swirl angle profile is input (angle in radians). 
Number of points in the spanwise profile (int). 

If > 0, dimensional spanwise positions used (feet). 

If < 0, fractional spanwise positions used (hub= 0.0 to shroud = 1.0). 
Inlet total pressure (real). 

Inlet total temperature (real). 

Inlet absolute flow angle (real). 

Spanwise position of inlet condition data (real). 
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4.5 Output File Description 


ADPAC-APES produces one or more output files depending on the parameters 
set in the input file. The user is referred to the previous sections for definitions of 
all the possible output that can be written. Of particular importance to restarting 
a solution or running a multiple blade row solution are the files listed in the Input 
File Description-Blade Row File Name Section. A flowfield restart, blade body force, 
and axisymmetric flowfield are written out using the Scientific DataBase Librarj 
routines [2]. These files are in a machine independent format, and are compatible 
with the PL0T3D multiple-grid, binary file description. 

Much of the output is written to the standard output device, which typically 
would be the system monitor, or can be redirected to a file in a unix system with the 
statement structure shown below. 

apes < input. file > output. file 

This statement will send the ADPAC-APES output to the file output. file. This output 
consists of a print out of the input parameters, messages of actions taken during 
the run (e.g., writing out a restart file), and a print out of the residuals and certain 
flowfield information at every iteration. An example of the “iteration history” portion 
of the output file is seen in Fig. 4.4. At the end of each iteration or step, some 
residuals of the flowfield variables are statistically manipulated and printed out. At 
each step the maximum value of the density residual (which can be thought of as a 
time derivative of the change in density) is printed out. The location of this maximum 
is also printed out and is sometimes useful in determining a region of instability in 
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TIME * 1279.3789 


SOLUTION RESTARTED FOR BLADE ROW 1 AT N - 300 


eddy ; 
step 

: cpu =2 . 1 

max (dr /dt) 

06373572 

i j 

seconds 
k rms(dr/dt) 

rssq(q) 

mf (in) 

mf (out) 

mf (core) 


301 

0 . 130E+02 

47 

54 

2 

0 . 251E+00 

0 . 279E+00 

289.73 

287.30 

38.32 


302 

0 . 126E+02 

47 

54 

2 

0 .252E+00 

0 . 280E+00 

289.74 

287.24 

38 . 32 


303 

0 . 123E+02 

47 

54 

2 

0 .251E+00 

0 . 279E+00 

289.74 

287 . 18 

38.31 


304 

0 . 119E+02 

47 

54 

2 

0 .251E+00 

0 . 278E+00 

289*. 7 4 

287 .11 

38 .31 


305 

0 . 115E + 02 

47 

54 

2 

0 . 251E+00 

0 .279E+00 

289.74 

287.05 

38 . 30 

i 

306 

0 . 110E+02 

47 

54 

2 

0 . 251E+00 

0.279E+00 

289.74 

286.98 

38 . 30 

i 

307 

0 . 105E+02 

47 

54 

2 

0 . 251E+00 

0 . 279E+00 

289.74 

286.92 

38.29 

| 

308 

0 . 100E + 02 

47 

54 

2 

0 . 252E+00 

0 . 279E+00 

289.74 

286.85 

38.29 

; 

309 

0 . 945E+01 

47 

54 

2 

0 . 253E+00 

0 . 280E+00 

289.74 

286.78 

38 .29 


310 

0. 891E+C1 

47 

54 

2 

0 .254E+00 

0 . 280E+00 

289.74 

286.72 

38.28 

= 


Figure 4.4: Portion of output file showing iteration history 


the flowfield or a poorly gridded region. The root mean square of the density residual 
is printed next. Following this the square of the change in all the state variables 
is summed and the root of this is printed (rssq(q)). The values printed after the 
residuals depend upon the geometry of the solution. If the geometry is unshrouded 
and there are no flow splitters, then a power coefficient, a thrust coefficient, and 
the number of supersonic flowfield points are printed. For shrouded geometries with 
no core defined, the inlet mass flow, exit mass flow, and the number of supersonic 
flowfield points are printed. For shrouded geometries with a core, the inlet, exit, and 
core corrected mass flow are printed (see Fig. 4.4). 

Three other output files are produced and they also have information on the 
residual history in them. They are plot files and a PostScript image can be generated 
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with the FULLPLOT program developed under the same NASA contract [4]. A 
description of FULLPLOT is reprinted from [4] in the Appendix. One file contains 
the history of the maximum density residual, another contains the history of the root 
mean square of the density residual, and the third contains both values. The program 
assigns a name to these files by adding a u ,his ( ” or u .mcizhisF or i .rmshist suffix to 
the fiowfield restart file name. For example, if the flowfield restart has the name 

flow. rest 

then the iteration history plot file of the density residual’s root mean square will be 
named 


flow. rest, rmshist 

Some of the output that the user can specify (through input parameters de- 
scribed in an earlier section) is printed in the nondimensional form used by the flow 
solver directly, and some output is printed out in dimensional form using the user set 
reference values. So that the user can understand and manipulate the output data 
better, a description of the flow variable nondimensionalization is presented. 

Some of the variables in the numerical solution are nondimensionalized by refer- 
ence values as follows: 


z r 

v z 

v r vq 

vr - , v e - 

v ref v ref 

T J ^ T 5 

L rtf L ref 

V Z — y 

v ref 

^ t 

II 

H 

* — K 

P 

P = „ > 

Pref 


The quantities are defined as follows: 
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p The nondimensional pressure. 
r The nondimensional radial length. 

T The nondimensional temperature. 

Vr The nondimensional radial velocity. 
vz The nondimensional axial velocity. 
vq The nondimensional tangential velocity. 
z The nondimensional axial length. 

p The nondimensional density. ; 

L re j The reference length nondimensionalizing the grid. \ 

This should be the input parameter DIAM. = 

p re j The reference (or freestream) relative total pressure 

This should be the input parameter PREF. | 

I 

T re j The reference (or freestream) temperature 

This should be the input parameter TREF. 

v Te j The reference (or freestream) velocity determined from the 1 

relative total conditions v re j = \JPrej!Pref 

R re f The reference (or freestream) gas constant 

This should be the input parameter RGAS. 

Pref reference (or freestream) relative total density 

determined from the relative total conditions 

Pref 

PTei = T re fR ref l9 
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4.6 Subroutine Description 


A list of the ADPAC-APES subroutines and their functions is given below for 
reference. A skeleton program flowchart is illustrated in Fig. 4.5. 

SUBROUTINE DESCRIPTION 


APES 

BAXI 

BCCORE 

BCEXT 

BCINL 

BCJV 

BCK 

BCKV 

BCSPL 

BCSPL2 

BCWAL2 

BCWALL 

BFORCl 


Main calling routine. 

Calculates axisymmetric flow quantities. 

Boundary condition for engine core simulation. 

Boundary condition routine for exit cells. 

Boundary condition routine for inlet cells. 

Boundary condition routine for j=constant surfaces (hub and 
outer boundaries). 

Inviscid boundary condition routine for the blade surfaces. 

Viscous boundary condition routine for the blade surfaces. 

Boundary condition routine for the hub, shroud, and flow splitter sur- 
faces. 

Axisymmetric boundary condition routine for the hub, shroud, and 
flow splitter surfaces. 

Axisymmetric boundary conditions at the exit for the 
2D solution. 

Boundary condition routine for periodic cells. 

Routine for calculating the axisymmetric average of the 
flowfield residuals. 
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BFORC2 

BFORC3 

BFORC4 

BFREAD 

CONATAN 

CONVAS 

EDDYDW 

EDYDW2 

ERROR 

EULER 

FILTER 

FIXUP 

FLREAD 

FSTR2 

FSTRES 

GETMAS 

GRIDG 

GSTR2 


Routine for calculating body force terms using the 
axisymmetric equations. 

Routine for calculating viscous terms needed in BF0RC1. 

Routine for calculating viscous terms needed in BF0RC2. 

Routine to read in body force terms for all adjacent blade 
rows. 

Continuous tangent function evaluation routine. 

Array storage conversion routine used in conjunction with SDBLIB 
machine-independent input/output routines. 

Turbulence model evaluation routine. 

Axisymmetric turbulence model evaluation routine. 
Convergence-checking routine for Runge-Kutta solver. 

Runge-Kutta solver. 

Artificial dissipation routine. 

Routine to correct corner phantom cells so output averages 
and viscous stress evaluations are correct. 

Routine to read in the flowfield restart file. 

Axisymmetric viscous stress evaluation along i coordinate direction. 
Viscous stress evaluation along i coordinate direction. 

Routine to compute mass average quantities at an i=constant coordi- 
nate surface. 

Routine to read and set up the grid. 

Axisymmetric viscous stress evaluation along j coordinate direction. 
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GSTRES 

HSTR2 

HSTRES 

INIT 

INITA 

INITB 

INITC 

INITD 

INPUT 

METRI2 

METRIC 

OUTPUT 

P3DOUT 

PLHIST 

PTCALC 

REFIN 

REFIN2 

REFIN2B 

REFIN3 


Viscous stress evaluation along j coordinate direction. 

Axisymmetric viscous stress evaluation along k coordinate direction. 
Viscous stress evaluation along k coordinate direction. 

Routine to initialize flowfield with reference flow quantities. 

Routine to initialize flowfield with axisymmetric solution. 

Routine to replace flowfield with a piecewise constructed one from the 
latest axisymmetric solutions. 

Routine to redefine flowfield from inlet to blade leading edge and from 
trailing edge to exit. 

Routine to reset inlet conditions based on the previous blade row ax- 
isymmetric solution. 

Routine to read in input data and set up reference values. 

Routine to calculate cell volumes and surface normals for 
2D solution. 

Routine to calculate cell volumes and surface normals. 

Routine to print output and save restart files. 

Routine to output PL0T3D flowfield (Q) files. 

Routine to output convergence history plot files. 

Routine to determine inlet boundary conditions. 

Routine to refine a 3D solution. 

Routine to refine a 2D solution. 

Routine to interpolate blade body forces to a finer grid. 

Routine to interpolate flowfield to a finer grid. 
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REFINB Routine to interpolate beta to a finer grid. 

RESID Implicit residual smoothing routine. 

RUNGE Runge-Kutta flux calculation routine. 

SAREA Routine to determine metric area terms. 

SAREA2 Routine to determine axisymmetric metric area terms. 
STEP Routine to determine time step. 

TRIM Routine to delete trailing blanks from a character string. 
STEP2 Routine to determine time step for 2D solution. 
VOLUX Routine to calculate an individual cell volume. 
VOLUX2 Routine to calculate an individual 2D cell volume. 

WF Wall function turbulent stress routine. 


It should also be mentioned that a number of routines from the SDBLIB library I 

i 

I 

are also included in the source code distribution, but are not defined here. 
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BF0RC1 


BCKV 


BCKV 


BF0RC2 


FIXUP 


PLHIST 


METRIC 


TRIM 


BF0RC4 


BCSPL2 


BCWAL2 


EDYDW2 


CONVAS 


FSTR2 

GSTR2 

HSTR2 


Figure 4.5: Subroutine flowchart of ADPAC-APES 
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5. ADPAC-APES SOLUTION TECHNIQUES 


5.1 Single/Multiple Blade Row Procedure 


A procedure for obtaining a numerical solution for multiple blade rows with the 
APES program is described below. The single blade row case is in general a reduction 
of the multiple blade row case and is described later. 

Before executing the solution algorithm, numerical grids (one for each blade 
row) are required. These grids model the actual three-dimensional geometry of their 
particular blade row, and represent the rest of the domain as an axisymmetric duct. 
APES requires all grids to have the same meridional representation (i.e., the same 
dimensional (z,r) coordinate lattice structure). If the geometry contains radial flow 
splitters (e.g., a cowl), the spanwise grid surface index is duplicated for the points 
defining the lower split surface. This is done so that the phantom cells defined in 
the flow solver will have a grid index associated with them. Figure 5.1 shows this 
spanwise index relationship. Grids produced by the TIGG3D program fulfill these 
requirements and more information on the grids and TIGG3D can be found in [5]. 

The solution procedure is begun from a set of initial data. This initial data is 
specified as a uniform flow (based upon the freestream Mach number), or may be 
introduced from a previous solution. The type of initialization selected, as well as all 
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J = 14 
J = 13 


J = 11, 12 


J = 10 


Figure 5.1: Spanwise index numbering around radial splitters 
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other program options, is controlled by parameters in the input file (see the Input 
File Description sections). For uniform flow initialization, the parameter REST must 
be zero and the parameter FSET should be zero as well. If a previous solution is 
used, the parameter REST is set to one. This procedure is typical for either a single 
blade row solution or the first blade row of a multiple row solution. For initialization 
of the blade row solutions following the first, the parameter REST is set to zero in 
those blade row input files. This will cause the initial flowfield to be computed from 
the axisymmetric solution of the previous blade row. Also, blade body forces from 
the previously run solutions will be incorporated in the current blade row calculation. 
The axisymmetric solution and body forces are read from file names listed in the input 
file. A diagram of this procedure for a three blade row solution is seen in Fig. 5.2. 

Once the initialization has occurred, a multiple blade row solution is found using 
a nested iteration procedure with an inner and outer loop as seen in Fig. 5.3 . The 
inner loop represents the iterative technique used to solve the flow equations for a 
particular blade row (a single APES run). This APES run for a individual blade row 
incorporates the effects of the neighboring blade rows (input through the body force 
file names). The run iterates a number of steps (marches in time to a steady solution) 
based on the parameter FNCMAX. For the APES restart/run the parameter REST 
must be set to negative one. When each blade row has gone through this iteration 
and the blade row effects have been recalculated, one cycle or “flip” through the 
system (outer loop) is complete. When these updated effects (which are calculated 
automatically) are used depends upon user control of the solution procedure. This 
control is achieved by what body force file names are put in the input file. For the 
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BLADE ROW #1 RUN 

FNCMAX = 100 



Freestream initialization 
(REST=0.) 

FILES READ 
None 


FILES WRITTEN 

flow.new.l bforc.new.l axi.new.l 


Restart from old solution 
(REST=1.) 

FILES READ 

flow.new.l bforc.new.l axi.new.l 


FILES WRITTEN 

flow.new.l bforc.new.l axi.new.l 


Initialize from previous 
blade row axi. solution 
(REST=1.) 

FILES READ 

flow.new.l bforc.new.l axi.new.l 

FILES WRITTEN 

flow'.new.2 bforc.new.2 axi.new.2 


Initialize from previous 
blade row axi. solution 
(REST=1.) 

HLES^READ 

flow.new.l bforc.new.l axi.new.l 
flow.new.2 bforc.new.2 axi.new.2 

FILES WRITTEN 

flow.new.3 bforc.new.3 axi.new.3 


Figure 5.2: Diagram of a three blade row solution initialization 
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example of an N blade row calculation (see Fig. 5.3), the first blade row’s input file 


could have a file name section as seen below. 


flow. old. 1 

bforc .old. 1 

axi . old. 1 

grid.l 

flow. old. 2 

bforc .old. 2 

axi .old. 2 

grid. 2 

flow. old. N 

bforc. old. N 

axi .old.N 

grid.N 

flow. new . 1 

bforc .new. 1 

axi .new. 1 



The bforc files on the first N lines are read in and used in the calculation. The 
bfovc name in the last line is the file where the updated blade effects for the first row 
are written. For the second blade row input file, the first line of the file name section 
could look like: 

flow. old. 1 bforc. old. 1 axi.old.l grid.l 

which would direct APES to use the first blade row’s body forces from the pre- 
vious flip. To use the updated blade row effects, the first line should instead look 
like: 

flow. new. 1 bforc. new. 1 axi.new.l grid.l 

For Fig. 5.3 the neighboring blade row effects are all from the previous flip, and 
the file names in the figure reflect that. However, if the updated terms for the current 
flip are used, the solution procedure is represented in 5.4. 

This nested loop procedure in continued until the convergence criteria is met. 
Typically, a solution is deemed converged when the root mean square average of the 
density residual has been reduced by a factor of 10 ^ from its initial value. For 
multiple blade row solutions, each flip is usually marked by an initial increase in the 
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Figure 5.3: Multiple blade row solution procedure 
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BLADE 

- ROW 
#1 


j 

\ 


BLADE 
- ROW 
#2 


J 

> BLADE 
ROWS 
3-N 


FILES READ 

flow.old.l bforc.old.l axi.old.l 
flow. old. 2 bforc.old.2 axi.old.2 

flow.old.N bforc.old.N axi.old-N 

FILES WRITTEN 

flow. new. 1 bforc.new.l axi.new.l 


FILES READ 

flow.new.l bforc.new\l axi.new.l 
flow.old.2 bforc.old.2 axi.old.2 

flow.oId.N bforc.old.N axi.oldN 

FILES WRITTEN 

flow.new.2 bforc.new.2 axi.new.2 


Figure 5.4: Solution procedure using updated blade row effects 
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residual, then a decrease so that a net reduction in the residual occurs. A second 
criterion (or measure of the trust worthness of the solution(s)) is that the mass flow 
at the inlet is close to the value at the exit. Typically, the difference between mass 
flows should be less than one percent of the inlet value. For multiple blade row 
solutions this criterion also applies to differences in mass flow between different blade 
row solutions. 

There is no well defined requirement for the number of flips needed before a 
solution is deemed converged. Many factors contribute (e.g., number of blade rows, 
strength of aerodynamic coupling) to the number of flips necessary for convergence. 
At least two flips is recommended as a minimum number and there is no recom- 
mended maximum. A general guideline is when the individual solution s residuals 
have dropped to an acceptable level, run one more flip. If the mass flow, pressure 
ratio, or RMS average of the residuals do not change significantly, the solution system 
should be converged. 

There are also no absolute rules governing the number of iterations each blade 
row must go through between flips. Generally, the number if iterations scales with 
grid size, with a reasonable choice being on the order of the number of axjal grid 
points in the solution. This number should be enough to let all primary perturbations 
run the coarse of the domain. A number much less than this might not allow this 
“sonic communication” to occur properly. Also a number much more than this can 
unnecessarily settle a solution before the whole solution system (all other rows) is 
at a steady state. During a solution startup (covered in the next section) or other 
significant transients, a large number of iterations on each blade row can sometimes 
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bring about unrealistic and unstable aerodynamic situations. 

For a single blade row case, the solution procedure is greatly simplified. With 
no neighboring blade rows and their effects to calculate, there is no outer loop. The 
iteration for the flow equation solver is executed until the convergence criteria is met. 


5.2 Solution Startup Techniques 

A solution startup procedure has been described for single and multiple blade 
rows in the previous section. The mechanics of the flow solver ( APES) were covered 
but nothing was mentioned about overcomming some of the difficulties with the fluid 
mechanics of a solution startup. This section suggests startup techniques first for 
single blade row, and then multiple blade row geometries. 

For a single blade row geometry, often a performance or flow condition is trying 
to be simulated. A match of rotational speed, mass flow, freestream Mach number, 
total pressure ratio, etc. is required. Setting the boundary conditions and other 
parameters to their final desired values and starting from freestream conditions might 
yield a converged, appropriate solution. However, with transonic fan sections, this 
often does not work. Starting up from freestream to an intermediate aerodynamic 
condition often yields good results. Two techniques are suggested and outlined below. 

Often getting the shock system “swallowed” in the supersonic regions of a blade 
row (see Fig. 5.5) is desired. One way to achieve this is to set the exit static pressure to 
a low value (below the inlet total pressure value) as an initial condition. This induces 
the flow to go in the right direction and results in the flowfield being in a choked 
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(maximum mass flow) condition. From that point, the exit pressure is increased to 
the final desired value. For the example of an isolated shrouded rotor, this procedure 
might look like path A in Fig. 5.6 on a performance map. 

A second procedure for starting up a transonic blade row is for the intermediate 
condition to be at a lower rotational speed than the final desired value. The rotational 
speed can be set to a fraction of the final value (.5 to .9) with the exit pressure set to a 
value appropriate to that condition. Once the solution is stabilized at this condition, 
the exit pressure and rotational speed are changed to their final values. This procedure 
is represented as path B on the performance map of Fig. 5.6. This procedure could 
eliminate the strong shock waves that occur with a full speed startup. These waves 
sometimes produce undesirable flow conditions and solution instabilities. 

For multiple blade row solutions, the same procedures for a single blade row are 
used. An intermediate condition is often desirable, but sometimes flow conditions re- 
quire a modified approach. For example, it is sometimes helpful to change the solution 
initialization process. The multiple blade row initialization procedure outlined in the 
previous section showed a continuous, serial progression through the rows. Another 
approach which has been successful is to initially make the outer loop iteration (see 
Fig. 5.3) between the first two blade rows alone. Once the two blade row system is 
stabilized, then the next blade row is initialized and included in the outer loop. This 
is repeated until all the blade rows are in the calculation. 
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Shock Unswallowed Shock Swallowed 


Figure 5.5: Shock system in a supersonic blade section 
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MASS FLOW 


Figure 5.6: Generic performance map for an isolated rotor 
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APPENDIX A. FULL PLOT : PostScript X-Y Plotting Routine 


The FULLPLOT program was developed to provide an automated plotting fa- 
cility for simple x-y type plots. The ADPAC-APES program produces files which 
monitor convergence history. These files are written in a format which FULLPLOT 
will use to create PostScript based plots. The execution of this program is extremely 
simple and can be invoked by the command: 


fullplot <filename 

where filename a plot file from an APES run. 

Array limits in the FULLPLOT program are determined by the PARAMETER 
statement: 


PARAMETER( MAXCUR=50, MAXPTS=5000) 


The parameter MAXCUR determines the maximum number of curves that FULLPLOT 
can plot, while MAXPTS determines the maximum number of points per curve. 

Input to FULLPLOT is directed upon execution, and this program can therefore 
easily be used for other two-dimensional plotting purposes. FULLPLOT produces a 
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PostScript plot file labelled fort. 15, which can be sent directly to a PostScript printer, 
or previewed on a compatible device. The FULLPLOT input file may be modified 
by the user to construct the final plots in a number of forms. A sample FULLPLOT 
input file is given in Fig. A.l. The PostScript plot for this input file is given in 
Fig. A.2. Each of the various line and symbol types, as well as the FULLPLOT 
grayscale shading are illustrated in the sample plot. 

The actual variables in the input file are free format. A description of each of 
the input PARAMETERS is given in the paragraphs below. 

VARIABLE DESCRIPTION 

TITLE Three lines of 80 character title data follow the TITLE header line. 

The 3 strings will be centered above the plot, and are plotted in the 
same order as they are input. Each title string must be in quotes. 

LVH Vertical/horizontal plot trigger: 

if = 1, the plot y axis is vertical on a standard 8^ x 11 page (portrait 
mode), 

if = 2, the plot y axis is horizontal on a standard 8^ x 11 page (land- 
scape mode). 

LEGTYP Legend type plot trigger: 

if = 0, no legend is plotted, 

if = 1, a plain legend is plotted, 

if = 2, a shadow box legend is plotted, 
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3 LINES MAXIMUM 

is the first title line' 
is the second title line' 
is the third title line' 

_ VERTICAL (0) OR HORIZONTAL (1) PLOT 


DXLEGI DYLEGI DXBCXI-- (DXLEGI , 

DXLOGOI DYLOGOI (LOGO -0 NO LOGO, DX, DY IN INCHES)* 

-1.3 -1-5 


T — . 

' This 
' This 
'This 
LVH — 

0 

LEGTYP 
2 

LOGC-- 
1 

X LABEL 

9 y a y ^ g ^ 

N T NCX- XSMIN XSMAX LOCXAX NXSIG- 

5 * 0.0 1.0 0 2 

YLA3EL 

'Y Axis : 

NINCY 

5 

NUMCUR— 

14 

CURVE da: 


DYLEGI INCHES FROM ORIGIN) 


,abel' 

-YSMIN 

-1 . 0 


— YSMAX LOCYAX NYSIG- 

3.0 0 2 


0.1 

C.l 

0.2 

C .2 

0.3 
C .3 

C.4 
C . 4 

C . 5 
0.5 

0.6 

0.6 

0.7 

0.7 

0.1 

0.1 

0.2 

0.2 

0.3 

0.3 

0.4 

0.4 

0.5 

0.5 

0.6 

0.6 

0.8 
C . 9 

o.e 

0.7 

0.8 


. ^ _ _ #PTS LT YPE , STYPE , LEGEND , GRAYP LT , GRAYSCALE , LEGLABE j. 

2 11 1*0 0.0 'Line 1 Symbol 1' 

0.1 
0 . 9 

2 2 2 1 0 0.0 'Line 2 Symbol 2' 

0 . 1 
0.9 

2 3 3 1 0 0.0 'Line 3 Symbol 3' 

0.1 

0.9 

2 4 4 1 0 0.0 'Line 4 Symbol 4 
0.1 
0.9 

2 5 5 1 0 0.0 'Line 5 Symbol 5' 

0 .1 
0.9 

26610 0.0 'Line 6 Symbol 6' 

0 .1 
0.9 

2 7 7 1 0 0.0 'Line 7 Symbol 7' 

0.1 

0.9 

2 6 B 1 0 0.0 'Line B Symbol 8' 

- 0.1 

- 0.9 

2 9 9 1 0 0.0 'Line 9 Symbol 9' 

- 0.1 

-0.9 

2 10 10 1 0 0.0 'Line 10 Symbol 10' 

- 0.1 

-0.9 

2 11 11 1 0 0.0 'Line 11 Symbol 11' 

- 0.1 

-0 . 9 

2 12 12 1 0 0.0 'Line 12 Symbol 12' 

- 0.1 

-0.9 

2 13 13 1 0 0.0 'Line 13 Symbol 13' 

- 0.1 

-0.9 


- 0.1 

-0.5 

-0.9 

-0.5 

- 0.1 


5 1 0 1 1 0.5 'Gray Shaded - Grayscale-0.5' 


Figure A.l: Sample input data file for FULLPLOT plotting program 
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Y Axis Label 


This is the first title line 
This is the second title line 
This is the third title line 


3.00 — | 


2.20 H 


i.40 H 



Line 1 Symbol 1 

Line 2 Symbol 2 

Line 3 Symbol 3 

Line 4 Symbol 4 

Line 5 Symbol 5 

Line 6 Symbol 6 

Line 7 Symbol 7 

Line 8 Symbol 8 

Line 9 Symbol 9 

Line 10 Symbol 10 

Line 1 1 Symbol 1 1 

Line 12 Symbol 12 

Line 13 Symbol 13 

Gray Shaded - Grayscale=0.5 



Figure A. 2: Sample PostScript x-y plot from FULLPLOT plotting program 
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DXLEGI 

DYLEGI 

DXBOXI 

LOGO 

DXLOGOI 

DYLOGOI 

XLABEL 

NINCX 


Legend placement variable. This variable determines the actual phys- 
ical length, in inches, from the lower left hand legend corner to the x 
axis. 

Legend placement variable. This variable determines the actual phys- 
ical length, in inches, from the lower left hand legend corner to the y 
axis. 

Legend box length parameter. The legend height is predetermined, 
but the legend box length is not. This variable permits user control of 
this function. This variable controls the actual physical length of the 
legend box in inches. 

Logo plotting trigger: 

if = 0, no logo is plotted, 

if = 1, the NASA logo is plotted. 

Logo placement variable. This variable determines the actual physical 
length, in inches, from the lower left hand logo corner to the x axis. 
Logo placement variable. This variable determines the actual physical 
length, in inches, from the lower left hand logo corner to the y axis. 

X axis label. This is an 80 character string which will be centered 
below the x axis. The string must be in quotes. 

Number of scale increments along the x axis. The x scale is determined 
by the number of increments and the minimum and maximum x axis 
scale values XSMIN, and XSMAX. The scale increment is then simply 
(XSMAX-XSMIN) / (NINCX-1) 
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XSMIN X axis minimum scale value (see NINCX). 

XSMAX X axis maximum scale value (see NINCX). 

LOCXAX Not used. 

NXSIG Number of significant digits past the decimal point in x axis scale 
markers. 

YLABEL Y axis label. This is an 80 character string which will be centered 
along the y axis. The string must be in quotes. 

NINCY Number of scale increments along the y axis. The y scale is determined 
by the number of increments and the minimum and maximum y axis 
scale values YSMIN, and YSMAX. The scale increment is then simply 
(YSMAX-YSMIN) / (NINCY-1) 

YSMIN X axis minimum scale value (see NINCX). 

YSMAX X axis maximum scale value (see NINCX). 

LOCYAX Not used. 

NYSIG Number of significant digits past the decimal point in y axis scale 
markers, 

NUMCUR Number of separate sets of data (curves) to be drawn. Each curve is 
represented by a set of x-y data points which can be represented by 
points, lines, or as a shaded region on the plot. 

The remainder of the FULLPLOT data set is the actual data to be plotted. The 
following data information must be repeated for each separate curve to be plotted. 
Each set of data (curve) begins with a header line with the following information: 

#PTS Number of x-y data pairs in the curve. 
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LTYPE Line type used to draw the curve. A curve may be represented 
by lines, symbols, shaded regions, or all of the above. 

The line types are defined as: 

LTYPE=0 is no line (symbol only), 

LTYPE=1 is a solid line, 

LTYPE=2 is a medium dashed line, 

LTYPE=3 is a short dashed line, 

LTYPE=4 is a medium chained dash line, 

LTYPE=5 is a very long dashed line. 

LTYPE=6 is a very short dashed line. 

LTYPE=7 is a staggered chain dashed line. 

LTYPE=8 is a alternate dash length line. 

LTYPE=9 is a very wide spaced dashed line. 

LTYPE=10 is a very short spaced short dashed line. 

LTYPE=11 is a tight staggered chained dashed line. 

LTYPE=12 is a very tight staggered chained dashed line. 
LTYPE=13 is a very short dash chained line. 

STYPE Symbol type used to draw the curve. A curve may be represented 
by lines, symbols, shaded regions, or all of the above. 

The symbol types are defined as: 

STYPE=0 is no symbol (line only), 

STYPE=1 is an open square, 

STYPE=2 is an open triangle pointing up, 
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STYPE=3 is an open triangle pointing down, 

STYPE=4 is an open circle, 

STYPE=5 is an black filled square, 

STYPE=6 is an black filled triangle pointing up, 

STYPE=7 is an black filled triangle pointing down, 

STYPE=8 is an black filled circle, 

STYPE=9 is an gray filled square, 

STYPE=10 is an gray filled triangle pointing up, 

STYPE=11 is an gray filled triangle pointing down, 

STYPE=12 is an gray filled circle. 

STYPE=13 is a horizontal bar. 

LEGEND A trigger to indicate whether the current curve should be included in 
the legend: if = 0, don’t include this curve in the legend, if = 1, include 
this curve in the legend. 

GRAYPLT A trigger to indicate whether the current curve should be plotted as 
a shaded region. If this is desired, the data will be interpreted as a 
closed curve which will have a shaded interior based on the shading 
value GRASCALE: 

if = 0, don’t treat this curve as a shaded region, 
if = 1, treat this curve as a shaded region. 

GRASCALE Shaded region shading value. This value must lie between 0.0 and 1.0. 

The value 0.0 represents true black, while a value of 1.0 represents true 
white. All other values in between represent increasing lighter shades 
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of gray. This value is only used when GRAYPLT— 1. 

LEGLABEL A character string (must be in quotes) indicating the legend label for 
this curve. 

The remaining lines represent the #PTS pairs of x,y data in free format. One 
data set is given per line for the given curve. 
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